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The critical densities and impact of forming A(1232) resonances in neutron stars are investigated 
within an extended nonlinear relativistic mean-field (RMF) model. The critical densities for the 
formation of four different charge states of A(1232) are found to depend differently on the separate 
kinetic and potential parts of nuclear symmetry energy, the first example of a microphysical property 
of neutron stars to do so. Moreover, they are sensitive to the in-medium Delta mass mA and the 
completely unknown A -p coupling strength g p a- In the universal baryon-meson coupling scheme 
where the respective A-meson and nucleon-meson coupling constants are assumed to be the same, 
the critical density for the first A _ (1232) to appear is found to be pfffl = (2.08±0.02)po using RMF 
model parameters consistent with current constraints on all seven macroscopic parameters usually 
used to characterize the equation of state (EoS) of isospin-asymmetric nuclear matter (ANM) at 
saturation density po. Moreover, the composition and the mass-radius relation of neutron stars are 
found to depend significantly on the values of the g p a and m/\. 

PACS numbers: 21.65.Ef, 21.65.Cd, 21.65.Mn, 26.60.Kp 


I. INTRODUCTION 

Understanding properties of A(1232) resonances in 
connection with possible pion condensation 0-1 in neu¬ 
tron stars and density isomers in dense nuclear mat¬ 
ter is a longstanding challenge of nuclear many- 
body physics. In fact, the role of A(1232) resonances 
in neutron stars has long been regarded as an impor¬ 
tant, and unresolved, issue j9j. Significant works have 
been carried out to understand in-medium properties of 
A(1232) resonances as well as their effects on satura¬ 
tion properties of nuclear matter and the equation of 
state (EoS) of dense matter using various many-body 
theories and interactions: see, e.g., refs. IToHTjl . How¬ 
ever, compared to the numerous investigations on the 
possible appearance and effects of other particles, such 
as hyperons and deconfined quarks, much less effort has 
been devoted to the study of A(1232) resonances in neu¬ 
tron stars in recent years. This is probably partially 
because of the rather high A(1232) formation density 
pX- i n the core of neutron stars pre dicted in the semi¬ 
nal work by Glendenning et al. [16l - |l8| using a mean-field 
model with parameters well constrained by the exper¬ 
imental data available at the time. Using default pa¬ 
rameters of their model Lagrangian leading to a sym¬ 
metry energy of E sym (p 0 ) = 36.8MeV and its density 
slope L(p 0 ) = 3pod E sy rn (p) /dp\ n=nn > 90MeV at satura¬ 
tion density po |16h 18| |. and using the universal baryon- 
meson coupling scheme in which the nucleon-meson cou¬ 
plings are set equal to the A(1232) -meson couplings 
(s<ta/3<tN = = ffpA/ffpN = 1), the critical den¬ 

sity p^‘4 above which the first A~(1232) appears is above 
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9po- This led to the conclusion that A(1232) resonances 
played little role in the structure and composition of neu¬ 
tron stars. In the same studies, the extreme importance 
of symmetry energy for the formation of both hyperons 
and A(1232) resonances was emphasized. Particularly, 
turning off the A -p coupling which contributes the poten¬ 
tial part of the symmetry energy (retaining thus kinetic 
symmetry energy only), the p)£if is about 3po [Tuj j. 

Interest has been renewed with recent studies using 
different symmetry energies and/or assumptions about 
the baryon-meson coupling constants which have found 
that the p™4 can be as low as po and the inclusion of 
the A(1232) has significant effects on both the composi¬ 
tion and structure of neutron stars mil. These studies 
generally use some individual sets of model parameters 
leading to macroscopic properties of ANM at saturation 
density consistent with most if not all of the existing ex¬ 
perimental constraints. 

During the last three decades, much progress has been 
made in constraining the EoS of dense neutron-rich nu¬ 
clear matter. In particular, reasonably tight constraints 
on the density dependence of the nuclear symmetry en¬ 
ergy E sym (p) especially around the saturation density 
have been obtained in recent years, see, e.g., refs. 1281431 
for comprehensive reviews. For example, the 2013 global 
averages of the magnitude and slope of the E sym (p) at 
po are respectively U sym (po) = 31.6 ± 2.7MeV and 
L = 58.9 ± 16.5 MeV based on 28 analyses of various ter¬ 
restrial laboratory experiments and astrophysical obser¬ 
vations [353 . Moreover, A(1232) resonances play a very 
important role in heavy ion collisions, see e.g., ref. [38} 
for reviews, especially for the production of particles 
such as pions, kaons and various exotic heavy mesons. 
In particular, the masses of A(1232) resonances primar¬ 
ily created in nucleon-nucleon (NN) collisions through 
the NN —> NA process act as an energy reservoir for 
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sub-threshold particle production. The release of this 
energy in subsequent collisions involving A(1232) reso¬ 
nances may help create new particles that can not be 
produced otherwise in the direct, first-chance NN colli¬ 
sions. Thus, particle production has been widely used 
in probing in-medium properties of A(1232) resonances. 
Since A(1232) resonances and nucleons have an isospin 
3/2 and 1/2 respectively, the total isospin is 1 or 2 for the 
NA while it is 1 or 0 for the NN state. Because of the 
isospin conservation, the A(1232) production can only 
happen in the total isospin 1 NN channel. Therefore, 
the abundances and properties of A(1232) resonances 
are sensitive to the isospin asymmetry of the system 
as neutron-neutron pairs always have an isospin 1 while 
neutron-proton pairs can have an isospin 1 or 0. 

Naturally, both heavy-ion collisions and neutron stars 
are places where the isovector properties and interactions 
of A(1232) resonances are expected to play a significant 
role. Indeed, useful information about the symmetry en¬ 
ergy of dense neutron-rich matter has been extracted 
from studyin g p ion and kaon productions in heavy-ion 
collisions [3l|, [32| . It is especially worth noting that the 
isovector (symmetry) potential of A(1232) resonances 
was recently found to affect appreciably the ratio of 
charged pions in transport model simulations of heavy- 
ion collisions at intermediate energies [H]. However, to 
our best knowledge, no quantitative information about 
the isovector interaction of A(1232) resonances has been 
extracted yet from any terrestrial experiments. On the 
other hand, there are strong indications from both the¬ 
oretical calculations and phenomenological model analy¬ 
ses of electron-nucleus, photoabsorption and pion-nucleus 
scattering that the A(1232) isoscalar potential Va (real 
part of its isoscalar self-energy Es) is in the range of 
—30 MeV + Vn < Va < Vn with respect to the nucleon 
isoscalar potential Vn [H}. The in-medium masses and 
widths of A(1232) resonances are also the focuses of many 
experimental and theoretical studies using various reac¬ 
tions and techniques. To our best knowledge, however, 
no clear consensus has been reached yet. For instance, 
from analyzing the (p, n) invariant masses in the final 
state of heavy-ion collisions at SIS/GSI energies, indica¬ 
tions were found for an approximately —60 MeV mass 
shift for A(1232) resonances at the freeze-out of about 
1/3 the saturation density [25}. However, photo absorp¬ 
tion data and some advanced model calculations found 
no evidence of significant in-medium A(1232) mass shift 
[HI 1 2l\ . It is thus exciting that new proposals to ex¬ 
perimentally study at FAIR/GSI the A(1232) resonance 
spectroscopy and interactions in neutron-rich matter are 
being considered by the NUSTAR Collaboration [40]. 

Similar to the appearance of any other new hadron 
above its production threshold in neutron stars, the ad¬ 
dition of A(1232) resonances will soften the EoS and in¬ 
fluence the composition of neutron stars [7(11 [m! . Because 
of charge neutrality, depending on the individual popu¬ 
lations of the four different charge states of A(1232) res¬ 
onances, the density dependence of the proton fraction 


in neutron stars may be modified. Then different cooling 
mechanisms sensitive to the proton fraction may come 
into play above certain densities. Moreover, the forma¬ 
tion of A(1232) resonances may also push up critical den¬ 
sities for the appearance of various hyperons [H| . As no¬ 
ticed earlier in the literature and emphasized in ref. [9[, 
there are many interesting questions regarding properties 
of A(1232) resonances in dense matter and their impact 
on observables of neutron stars. Obviously, answers to all 
of these questions naturally rely on the critical density of 
A(1232) formation in dense neutron star matter. 

In this work, we first identify analytically key micro¬ 
physics quantities determining the critical formation den¬ 
sities of the four charge states of A(1232) resonances. 
Then, within a nonlinear RMF model we calculate con¬ 
sistently the p^'J as a function of the A(1232) mass toa, 
the isovector A-p coupling strength g p a and seven macro¬ 
scopic variables characterizing the EoS of ANM at po all 
within their latest constraints. Finally, impacts of the 
A(1232) formation on the composition and mass-radius 
correlation of neutron stars are studied. 


II. KEY MICROPHYSICS DETERMINING THE 
DELTA FORMATION DENSITY IN NEUTRON 
STARS 



P (fm‘ 3 ) 


Figure 1: The maximum mass of Delta resonances produced 
in the NN — > NA process in a free Fermi gas of nucleons at 
density p. 

To set a reference for our following studies we first es¬ 
timate the p^ 1 / in a free Fermi gas of nucleons. For the 
head-on collision of two nucleons both with Fermi mo¬ 
mentum |k| = fcp = (37r 2 p/2) 1 / 3 in the NN — > NA pro¬ 
cess, the maximum mass of the produced A(1232) reso¬ 
nance is m^ ax = 2(kp + m 2 j) 1 ' 2 — ton where ttin is the 
average nucleon mass in free-space. It is well known that 
A(1232) has a Breit-Wigner mass distribution around 
the centroid rn° A = 1232 MeV with a width of about 120 
MeV. The distribution starts at a minimum of m™ n = 
TonJ r ~ 1076 MeV where m w is the pion mass. Shown 
in Fig. |T] is the m^ ax reachable in the NN —>- NA process 
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as a function of density. From this simple estimate, where 
effects of the nuclear potentials are neglected, the critical 
density p™ 1 for producing the lightest A(1232) resonance 
is about 3po. To reach the centroid = 1232 MeV, the 
density has to be far above lfm~ 3 . This estimate also 
illustrates the importance of considering the mass depen¬ 
dence of the A(1232) formation density in more realistic 
calculations for matter in neutron stars. 

In interacting nuclear systems, the masses of 
A(1232) resonances and their critical formation densities 
depend on the in-medium self-energies of all particles in¬ 
volved. Assuming neutron stars are made of neutrons, 
protons, A(1232) resonances, electrons and muons, i.e., 
the npepA matter in chemical and /3 equilibrium, the 
total baryon number density is 

p = Pp + Pn + Pa++ + Pa+ + Pa° + Pa - (1) 


where 


Pn = 


Pa++ = 


P A+ = 


(fcg) 3 (££) 3 

37T 2 ’ Pp 3 tt 2 ’ 

(^ ++ ) 3 „ _ MJ ! 

9 1 P A 9 

7T Z 7T Z 

(^ + ) 3 (kff 

3tt 2 ’ PA ° 3t r 2 


( 2 ) 

(3) 

(4) 


with fcp’s (j = p,n, A ++ , A + , A°, A - ) being the corre¬ 
sponding Fermi momenta. The chemical equilibrium con¬ 
dition for reactions n —»• p + e + V e and p + e —» n + v e 
requires p e = p n — p P where p e = [ml + (&!) 2 ] 1 ' 2 = 
[m 2 + (37r 2 pa;e) 2 / 3 ] 1 / 2 — (3n 2 px e ) 1 ^ 3 with x e = p e / p the 
electron fraction. When the chemical potential of elec¬ 
tron is larger than the static mass of a muon, reactions 
e— > p + i/ e +Vfj,, p + p ^ n + and n — > p + p + will 
also take place. The latter requires 


Pn Pp — Pfi — 



(37r 2 px M ) 2 / 3 


(5) 


besides p n — p p = p e , where = 105.7MeV is the 
mass of a muon and x^ = p^/p is the muon fraction. 
On the other hand, the following four types of inelastic 
reactions will take place between nucleons and the four 
charge states of Delta resonances 

A++ + n<—>p + p, (6) 

A++n<—>n + p, (7) 

A 0 + p <■—> p + n, (8) 

A - + p <—> n + n. (9) 


Their chemical equilibrium requires then 


Pa++ —2p p /i n , (10) 

Pa+ = Ppi ( 11 ) 

Pa° = Pm ( 12 ) 

Pa~ = 2 p n Pp- ( 13 ) 


In addition, the total charge neutrality in neutron stars 
requires that x p + x&+ + 2xa++ = x e + x M + x^- where 
xa~ = Pa~ / Pi = Pa + / Pi and xa++ = Pa++ / Pi 
respectively. 

Within the framework of a given nuclear many-body 
theory, the Eqs. m-m can be used to calculate the 
critical formation densities for the four charge states of 
A(1232) resonances. Generally speaking, in relativistic 
mean-field models, a baryon of bare mass Wbaryon obtains 
a Dirac effective mass mj irac (baryon) = TObaryon + ^s and 
a chemical potential pbaryon = [fcp + m))? rac (baryon)] 1 / 2 + 
Sy where Eg and Ey are the real parts of its scalar and 
vector self-energies, respectively. Consider the A - for¬ 
mation, for example, noticing that p n — p p ~ 4E synl (p)5 
and using non-relativistic kinematics, the Eq. m leads 
to the following condition for producing a A - of bare 
mass mA- at rest 


(A#) 2 


2 TO 


* 

dirac 


[37t 2 (1 + (5)p/2] 2/3 


2 TO 


* 

dirac 


(14) 


— m-A- ~ m N — 4E sym (p)<5 + E^ — Eg + Ey — Ey 


where m^ ilac I s the nucleon Dirac effective mass and 
8 = (p n — Pp)/p is the isospin asymmetry of nucleons be¬ 
fore A(1232) resonances are produced. Given the density 
dependencies of the symmetry energy and self-energies, 
this equation determines the p^_ in neutron stars at 
/3 equilibrium. It also shows clearly what microphysics 
quantities determine the p™! . In particular, the differ¬ 
ence in Delta and nucleon masses toa- — Ton, the symme¬ 
try energy £' sym (p) , the difference in both scalar E^ —Eg 
and vector Ey — Ey self-energies are all characteristics 
of baryon interactions. It also indicates where the model 
dependence and uncertainties are. As we mentioned ear¬ 
lier, experimental data indicates that the difference in 
nucleon and A(1232) isoscalar self-energies can be up to 
about 30 MeV while there is simply no experimental in¬ 
dication so far about the difference in their isovector self¬ 
energies. We notice that in many studies in the literature 
the A(1232) resonances and nucleons are assumed to have 
the same scalar and vector self-energies. In this case then, 
the p“‘i is completely determined by the A(1232) mass 
toa and the nuclear symmetry energy as a function of 
density E sym (p) . 

The nonlinear RMF model has been very successful 
in describing many nuclear properties and phenomena 
during the last few decades, see e.g., refs. |4ll - l59| . The 
total Lagrangian density of the nonlinear RMF model of 
ref. [54[ augmented by the Yukawa couplings of the Delta 
fields to various isoscalar and isovector meson fields, can 
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be written as M, m, 0 , mu 

£ =^N [ 7 ~ 9wNW M - SpNTN • p M ) 

- (ton - 5<rNCr)] 'i/j N 
+ t'W [7 m(*^ - S'^AW'" - ffpATA • p M ) 

- (TOa -ff CT ACr)]^A 

+ ^3 M cr9 M cr - ^to^ct 2 - {7 N (er) 

+ ^TO 2 ^^ - + ic^N^N ^'') 2 

.1 2 -* ->u ^ -» -*uz> 

+ 2 m p^M • ^ 4 ^ * P 

+ 2 (ffpNPp • p'") A V5 2 N W M W, "> ( 15 ) 

where = dpWj, - d v and p pl/ = d p p v - d v p p 
are strength tensors for w and p meson fields, respec¬ 
tively. i/'n, V’aj ct > w /ji Pp are the nucleon Dirac field, 
Schwinger-Rarita field for A resonances, isoscalar-scalar 
meson field, isoscalar-vector meson field and isovector- 
vector meson field, respectively, and the arrows denote 
isovectors, Un(<t) = b (T NmN(gaNcr) 3 /3 + c r7 N(g C rNcr) 4 /4 is 
the self interaction term of the a field. The parameter Ay 
represents the coupling constant of mixed interaction be¬ 


tween the isovector p and isoscalar w mesons. It is known 
to be important for calculating the density dependence 
of the symmetry energy [54| . In terms of the expecta¬ 
tion values of the meson fields, cf, uJo and Pq 3 * where 
the subscript “0” denotes the zeroth component of the 
four-vector while the superscript “(3)” denotes the third 
component of isospin, the nucleon and A(1232) isoscalar 
self-energies are respectively 

s S = -So-nct (16) 

and = —g a act. (17) 

Their isovector self-energies are respectively 

£y = SwNWo + T p/n5pNPo ^ (18) 

and £y = g u awo + T 3 g p APo ^ (19) 

with Tp = +1, r 3 = -1 and i = A++, A+, A 0 , A - , 

T 3 — -4-3 T 3 — -1-1 T 3 — —1 T 3 — —3 
'A++ — wo, / A+ — -|-i, / A o — ' A - — a. 

In terms of the ratios of Delta-meson over nucleon- 
meson coupling constants x a = g a A/gaN, x^ = P^a/p^n 
and x p = PpA/ffpN, the Eas. (flQl) - (fl3l) lead to the following 
conditions determining the critical densities for forming 
the four charge states of A(1232) resonances 


J 


„crit . {*%) 


PA- '■ 


Jm * 
Z ' U dirac 


=4$A + Po-n(1 - x a )a - 5wn(1 - Xu,)ujo - 6(1 - x p )E p °^(p)6 - 4.E k “ (p)c). 


Pao ■ $a + Po-n(1 - X„)a - 5ojn(1 - x w )ujo - 2(1 - x p )E p °^(p)S : 

^ TO dirac 
„crit . (*?) 


Pa+ : 


^ TO dirac 


^ $A + Qtj N (1 - X a )a - 5taN(l - Xuj)loo + 2(1 - Xp)E^(p)S 1 


(k p ) 

" crit • —~d>A+PrrN(l-a; I T)a : -5a;N(l-a4a;)wo-|-6(l-a;p)£; s p y 0 ^(p)5 + 4£; s k 3 )(( 1 (p) 

Zm dirac 


Pa++ : 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 


where 4 >a = toa — Ton is the Delta-nucleon mass difference, E p °^(p) = 2 1 pg 2 p ^{ , m 2 p + Avp 2 n p^ n Wq) 1 and E^ l n (p) = 
6 -1 fc 2 (A: 2 + to*; 2 „„)~ 1//2 are respectively the potential and kinetic part of the symmetry energy in the nonlinear RMF 
model |32|. 

I- 


Several interesting conclusions can be made qualitatively 
from inspecting the above four conditions. Generally, 
the critical densities depend differently on the three cou¬ 
pling ratios x a , x u and x p as they have different na¬ 
tures. The isoscalar coupling ratios x a and x u affect the 
four A(1232) resonances the same way, i.e., the x a low¬ 
ers while the x^ raises their critical formation densities, 
while the isovector coupling ratio x p acts differently on 
the four different charge states of A(1232) resonances. 
Moreover, the kinetic and potential parts of the symme¬ 
try energy have separate and different effects. In par¬ 
ticular, in the universal baryon-meson coupling scheme, 
i.e., x a = x^ = x p = 1, the critical densities for cre¬ 


ating A and A ++ depend only on the <I>a and the ki¬ 
netic symmetry energy (p) besides the m)j irac , while 
those for the A° and A + are determined only by the 
<!>a and mj irac . In this case, assuming the E^^p) is 
always positive as in the case of RMF, noticing that the 
Fermi surface of protons is lower than that of neutrons 
at any density in neutron-rich matter, i.e., kp < kp, one 
then sees immediately the following sequence of appear¬ 
ance p c £l < p A o* < p A i t < p A ) t + However, we 

notice that the short-range nucleon-nucleon correlation 
(SRC) [fHU may lead to negative kinetic symmetry en¬ 
ergies even at normal density of nuclear matter, see e.g. 
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refs. [63l469j . In this case, the order of appearance of A - 
and A ++ , thus the fraction of various particles and the 
structure of neutron stars may be different. We remark 
here that this is the first time that some physics quanti¬ 
ties in neutron stars are found to depend separately on 
the kinetic and potential parts instead of the total sym¬ 
metry energy. 



Figure 2: (Color Online) Dependence of the critical density 
p c £- for A“ formation in neutron stars on the relative A — p 
coupling strength x p = s^a/s^n- 

Some earlier studies, see, e.g., refs. @, @, [H, [24[ , indi¬ 
cate that x a « Xu « 1.0. Effects of slight deviations from 
this value on properties of neutron stars have also been 
reported, see, e.g., ref. Hi]. To our best knowledge, how¬ 
ever, little is known about the range of x p and its effects 
in either heavy-ion collisions [39,] or neutron stars. More¬ 
over, most studies so far are limited to density isomers 
due to A(1232) formation in symmetric nuclear matter 
where it is sufficient to consider only effects of the x a and 

Xjj. 

In the universal baryon-meson coupling scheme, the 
default set (SH-NJ) of RMF model parameters leads 
to the following values of macroscopic quantities char¬ 
acterizing the EoS of ANM at po = 0.149 fm -3 : the 
binding energy E 0 (p 0 ) = —16.09MeV, the Dirac effec¬ 
tive mass TOji rac (po)/?7iN = 0.64, the incompressibility 
Ko(po) = 230 MeV, the skewness coefficient Jo{po) = 
—415 MeV, the magnitude E sym (p 0 ) = 31.17MeV and 
slope L(pa) = 48.64 MeV of symmetry energy. We re¬ 
mark that almost all of these bulk parameters are ex¬ 
tracted from a fit to the properties of finite nuclei jUj 
with the exception of the skewness coefficient Jq that 
has been obtained by fixing the maximum mass of a neu¬ 
tron star at M max = 2.01 M 0 [70] . Variations around this 
parameterization will be investigated in the following. 

For Delta resonances of mass toa = 1232 MeV, we 
study in Fig.[2]the p™d dependence on the x p while keep¬ 
ing all other quantities at their default values. It is seen 
that the p^d increases approximately linearly with x p . 


At x p = 1, the p^d is only about 2.1po consistent with 
that found in refs. [2l|, [23, [24| but much smaller than the 
one found by Glendenning et al. However, un¬ 

less the value of x p is somehow constrained, the p^d will 
remain underdetermined. 



m A (MeV) 


Figure 3: (Color Online) The Delta mass toa dependence of 
the critical density p^d for A - formation in neutron stars 
(blue) and the Breit-Wigner mass distribution of Delta reso¬ 
nances in free-space (red). 

Delta resonances in free-space have the Breit-Wigner 
mass distribution 

/(mA) = 4 (toa — + r 2 (?7iA)/4 (24) 

with the mass-dependent width given by [nUzl 

r(m A ) = 0.47 9 3 /(u^ + 0 . 6 q 2 ) (GeV) (25) 

where q = [([m\ — + ?71^]/2toa) 2 — to 2 ] 1 / 2 is the 

pion momentum in the A rest frame in the A —> w + N 
decay process. Shown in red in Fig. [3] is the free- 
space A(1232) mass distribution. It is known that the 
A(1232) mass distributions may be modified in nuclear 
medium [Z(|. This effect is beyond the scope of the RMF 
model used here as it does not consider the imaginary 
part of the A(1232) self-energy self-consistently. How¬ 
ever, we can examine how the p^d depends on the bare 
A(1232) mass toa by varying its value in the Eqs. E0l) - 
G51) . As one expects and indicated by the Eqs. C0l) - 
(123[) . the p™d increases with toa in the universal cou¬ 
pling scheme. Our numerical calculations shown with 
the blue line indicate that the increase is almost linear. 
Considering the mass distribution, while the p^d for A 
resonances around m a is about 2.1po, it gradually de¬ 
creases for lower A masses. Of course, these low-mass 
A resonances are less likely to be produced compared to 
the ones near On the other hand, the A mean life¬ 
time ta = h/T(rn&) is only about 1.7 fm/c at but 
increases very quickly for lower masses. Thus, the main 
population of A resonances in neutron stars may not nec¬ 
essarily peak at toa = 1232 MeV. A detailed study of 
this issue will require a full account of the 7r — N — A 
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dynamics in neutron stars that is also beyond the scope 
of the RMF model used here. Nevertheless, our results 
indicate that the appearance of A resonances, especially 
the ones with low masses around 2po, may compete with 
other particles, such as hyperons, thus possibly modify 
the widely accepted and longtime viewpoint that hyper¬ 
ons should appear earlier than A(1232) resonances in 
neutron stars [3, [73}. To our best knowledge, however, 
no study to date has considered consistently effects of 
the mass distribution and the associated mass-dependent 
lifetimes of A resonances in neutron stars. 


III. EFFECTS OF NUCLEAR EQUATION OF 
STATE ON THE FORMATION OF 
A(1232) RESONANCES IN NEUTRON STARS 


The equation (fl3l) for determining the critical density 
can be rewritten as 


Ma- = to a + S a = 2/i n - p p ~ p n + 4£’ sym (p)<5. (26) 


The left hand side is given in terms of the micro¬ 
scopic quantities, i.e., A(1232) mass and the three Delta- 

meson coupling constants in the total self-energy E A = 
_ _(3) 

-ffo-AU + g u AWo - 3p p AP 0 ■ Using the parabolic ap¬ 
proximation for the EoS of ANM E(p,5 ) ~ Eq(p) + 
E syin (p)S' 2 +0(S 4 ) and assuming the density is not too far 
from p 0 , the right hand side of Eq. (l26l) can be expanded 
in terms of y = (p — po)/3po and the isospin asymmetry 
6 as 


Ma*- — to n + -E'o(po) + 


XP_ 
3 Po 


X_ 

2 


K 0 


X 


P 


+ X ) Jo + (25 + 3S 2 )E sym (po) 


6 \Po 

7r—d 2 + x(2<5 + 3<5 2 ) L(po) 
,3p 0 


(27) 


where higher order terms in the expansion have been 
neglected (for full details see ref. (58(). Here Kq = 
9pod 2 E 0 (p)/dp 2 | p=p? and J 0 = 27pgd 3 £ , o (70/dp 3 | p=po 
are the incompressibility and skewness coefficient of sym¬ 
metric nuclear matter (SNM) at po, respectively. This ex¬ 
pansion is very easy to understand considering the energy 
conservation in the Delta production process NN -A NA, 
i.e., the minimum energy of the Delta is the energy of 
a nucleon (sum of nucleon rest mass and its mechani¬ 
cal energy). Since all seven macroscopic quantities used 
to characterize the EoS of ANM i.e., (a) the satura¬ 
tion density p 0 of SNM where the pressure P(po) = 0, 
(b) the binding energy Eo(po), (c) incompressibility Kq, 
(d) skewness coefficient Jo, (e) nucleon effective mass 
TO d?rao ( f ) magnitude E sym (p 0 ) and (g) slope L(p 0 ) of 
the symmetry energy, are all explicit functions (see, e.g., 
refs. [53, |7J, |75[) of the seven RMF microscopic model 
parameters, i.e., g aN , g uN , g p n, b a N , c ctN , c wN and A v , 
the Eq. (ET71) allows us to explore the p™i dependence 


on either the seven microscopic or macroscopic param¬ 
eters. Since the macroscopic quantities are either em¬ 
pirical properties of nuclear matter or directly related to 
experimental observables, it is more useful for the pur¬ 
poses of this work to examine the p c £l by varying in¬ 
dividually the seven macroscopic quantities. We notice 
that within both the RMF and Skyrme-Hartree-Fock ap¬ 
proaches similar correlation analyses [76| have been suc¬ 
cessfully applied to study the neutron skin IE l77ll. th e gi¬ 
ant monopole resonances (GMR) of finite nuclei I78| . the 
higher-order bulk characteristic parameters of ANM [79| , 
the electric dipole polarizability qd in 208 Pb [§Gj], the cor¬ 
relation between the maximum mass of neutron stars and 
the skewness coefficient of the SNM j58j, as well as the re¬ 
lationship between the E sym (p) and the symmetry energy 
coefficient in the mass formula of nuclei [8l| . 

In Fig. [4] we show the p™d obtained in the universal 
coupling scheme as a function of the seven macroscopic 
parameters within their respective uncertain ranges. 
Among them, the Jq = —250 ± 250 MeV and L = 
50 ±20 MeV currently have the largest uncertainties. Ex¬ 
amining the results shown in Fig. 0 we make the follow¬ 
ing two observations: (i) The values of p™! from all seven 
correlations overlap around p™_( =(2.08 ± 0.02)po. No¬ 
tice that this result depends on the value of Jq in the 
default set, which had been fixed to obtain the maxi¬ 
mum stellar mass configuration matching the observed 
2.01Mq neutron star. With this assumption then this re¬ 
sult is the most reliable prediction for the p™d consistent 
with all of the constraints within the RMF model con¬ 
sidered. We note that increasing Jq raises the maximum 
mass and lowers the critical density p c £l ; hence the range 
p c £- =(2.08±0.02)p o constitutes an upper limit on p™i . 
Since the inclusion of Delta resonances inevitably softens 
the EOS, the magnitude of Jq should therefore be much 
larger than the one used in our default set to be consis¬ 
tent with the current observation of the two solar-mass 
neutron star. With the other six macroscopic quantities 
fixed within their current uncertainty ranges, then the 
critical density p™i will take an even smaller value than 
the upper limit predicted above. And this is one of the 
main reasons why we haven’t considered hyperons in the 
present work, because they appear at much higher densi¬ 
ties than p™i . (ii) A reasonable and quantitative mea¬ 
sure of the sensitivity of p c £l to each individual variable is 
the response function Q = \(dy / y) / {dx / x)\ where d y/y 
is the relative change in p with respect to its mean 
value and dx/x is the relative change in the variable x 
with respect to its mean value. The value of Q is approx¬ 
imately Q[Eq(pq)} = 0.53, <2[p 0 ] = 0.40, Q[K 0 ] = 0.29, 
2KLJ = 0T8, Q[J 0 ] = 0.17, Q[U sym (p 0 )] = 0.11 and 
Q[L(p 0 )\ = 0.10 following the same order of importance 
in the expansion of the minimum Delta chemical poten¬ 
tial in Eq. (fT7p . 

The seemingly stronger correlation between the p c £ If 
and Jq compared to the correlations of p c £} with the 
other 6 variables is because of the relatively larger un¬ 
certainty in Jq. On the other hand, some of the weaker 
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Figure 4: Critical density for the formation of A resonance from the nonlinear RMF model by varying individually po (a), 
Eo(Po) (b), m3Lc (c), K 0 (d), J 0 (e), £ sym (p 0 ) (f) and L(p 0 ) (g). 



L(wp 0 )/w (MeV) 

Figure 5: (Color Online) Critical density for the formation 
of A - resonance as a function of the scaled density slope 
of symmetry energy L(wpo)/w at three reduced densities of 
p r /po = w = 0.7, 1 and 2, respectively 


correlations shown in Fig. [4] may become much stronger 
if one goes beyond their current constraints. For exam¬ 
ple, shown in Fig. [5] are the correlations of the p^ 1 * with 
the reduced slope L(p r )/w of E sym (p) at three reference 
densities of p r /po = w = 0.7,1.0 and 2 while fixing the 
magnitudes of E sym (p) and other variables at their de¬ 
fault values. It is seen that the p™i increases much faster 
for L(p T )/w > 70 MeV. This feature is consistent with 
that found in ref. (24j. Noting again that some earlier 
studies have used models predicting L values high than 
90 MeV, they thus predicted correspondingly much large 
values for the p^l 1 . 


IV. EFFECTS OF DELTA RESONANCES ON 
THE COMPOSITION AND STRUCTURE OF 
NEUTRON STARS 

Having shown that the p^i! depends sensitively on the 
completely unknown A-p coupling strength g p a and the 
Delta mass toa, and it is about p™i =(2.08 ± 0.02)po in 
the universal coupling scheme for toa = 1232 MeV us¬ 
ing RMF model parameters consistent with all existing 
constraints on the nuclear EoS, we now turn to effects 
of Delta formation on properties of neutron stars. This 
study is carried out within the npepA model omitting 
other particles such as hyperons and quarks at high densi¬ 
ties. This model is sufficient for the purposes of this work. 
Moreover, we restrict ourselves to studying effects of the 
Delta formation on the composition and mass-radius re¬ 
lation of neutron stars by varying the Delta mass and 
its coupling strength with the p meson. In constructing 
the EoS of various layers in neutron stars for solving the 
Oppenheimer-Volkoff (OV) equation, we follow a rather 
standard scheme. For the core we use the EoS of /?- 
stable and charge neutral npepA matter obtained from 
the nonlinear RMF model described earlier. The inner 
crust with densities ranging between p out = 2.46 x 10 -4 
fm~ 3 corresponding to the neutron dripline and the core¬ 
crust transition density p t determined self-consistently 
using the thermodynamical method [dj. 82| is the region 
where complex and exotic nuclear structure —collectively 
referred to as the “nuclear pasta” may exist. Because 
of our poor knowledge about this region we adopt the 
polytropic EoSs parameterized in terms of the pressure 
P as a function of total energy density e according to 
P = a + 6e 4 / 3 0, [83| . The constants a and b are de¬ 
termined by the pressure and energy density at pt and 
Pout 0. For the outer crust 0, we use the BPS EoS 
for the region with 6.93 x 10 -13 fm -3 < p < p ou t and the 
FMT EoS for 4.73 x 10" 15 fm" 3 < p <6.93 x lO" 13 fm” 3 , 
respectively. 

Shown in Fig. [6] are fractions of different species, i.e., 
= Pi/Pi in neutron stars using two parameter sets 
with the A-p coupling strength corresponding to x p = 1 
















P ; Po 

Figure 6: (Color Online) Fractions of different species in neu¬ 
tron stars composed of neutrons, protons, Delta resonances, 
electrons and muons within the nonlinear RMF model using 
two sets of model parameters with x p = g p A/g P N = 1 and 2, 
respectively. 

and 2, respectively. These calculations are done with 
toa = 1232 MeV. It is seen that the appearance of A - 
affects significantly the fractions of others particles de¬ 
pending on the value of the x p as one expects. The mod¬ 
ified fractions of the lighter particles e and /i will affect 
their weak decays and thus the possible kaon condensa¬ 
tion. Moreover, the strong boost of the proton fraction 
above the A - production threshold may have a signifi¬ 
cant impact on the cooling processes in neutron stars [85j . 
A detailed investigation of these consequences requires a 
self-consistent extension of the model considered here and 
is on the agenda of our future work. 



Figure 7: (Color Online) The mass-radius correlation of neu¬ 
tron stars without and with Delta resonances of different 
masses using x p = <? p a/<?pN = 1 and 2, respectively. 

Finally, we study in Fig. [7] effects of both the Delta 
mass m a and the x p parameter on the mass-radius rela¬ 
tion of neutron stars. Without including the Delta reso¬ 


nances (black line), as we mentioned earlier the default 
value of Jo was chosen to predict a maximum neutron star 
mass of 2.01 Mq consistent with the latest observations. 
As given in Eq. EH) and Eq. E5l) . A(1232) has an in¬ 
trinsic mass distribution and each mass has its own width 
(lifetime) in free space. Especially, it is possible to form 
A(1232) with masses much smaller than the centroid of 
1232 MeV in low energy NN —> N A or 7 t + N —» A 
reactions. As one expects, including Delta resonances 
reduces both the maximum mass and the corresponding 
radius. Moreover, the effect is stronger for Delta reso¬ 
nances with masses smaller than the centroid. These low 
mass Delta resonances have lower production thresholds 
and are thus more abundant. They then softens the EoS 
more starting at lower densities reached in the cores of 
lighter neutron stars. Of course, effects of the Delta res¬ 
onates on properties of neutron stars depend on the value 
of x p . As shown in Fig. [G] and discussed earlier, a larger 
value of x p leads to a higher p c £If . Thus, the effect of 
varying toa appears at a higher mass (higher central den¬ 
sity) of neutron stars when the value of x p is increased 
from 1 to 2. 

Obviously, including the Delta resonances we can no 
longer obtain a maximum neutron star mass of 2.01 Mq 
with the default parameter set. Moreover, we found that 
by varying the value of Jo in its uncertain range, we can 
not recover a maximum mass of 2.01 Mq with x p = 1. We 
confirm that indeed there is a “Delta puzzle” consistent 
with the finding of refs. [23j, [24j . Thus, if indeed Delta 
resonances can be formed around 2po as in the universal 
coupling scheme, our findings call for serious considera¬ 
tions of the various EoS stiffening mechanisms proposed 
in the literature. On the other hand, as we mentioned 
earlier, there is currently no constraint on the x p neither 
theoretically nor experimentally. Since the critical den¬ 
sity for Delta formation increases approximately linearly 
with x p , it is therefore critical to first constrain indepen¬ 
dently the x p , such as from heavy-ion collisions, before 
the “Delta puzzle” in understanding properties of neu¬ 
tron stars can be resolved. 


V. SUMMARY AND DISCUSSIONS 

The possible formation and roles of A(1232) reso¬ 
nances in neutron stars are outstanding issues in nuclear 
astrophysics. The first and most important piece of in¬ 
formation necessary for resolving these issues is the crit¬ 
ical density above which the A(1232) can be formed in 
neutron stars. Previous studies have indicated that the 
critical densities range from po to very high values only 
reachable in the core of very massive neutron stars. In 
this work, within the extended nonlinear RMF model we 
found that the critical formation densities for the four 
different charge states of A(1232) resonances depend dif¬ 
ferently on the separate kinetic and potential parts of 
the nuclear symmetry energy, the in-medium Delta mass 
to a and the completely unknown A -p coupling strength 
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g p a. This is the first time a microphysical property of 
neutron star matter has been shown to depend differ¬ 
ently on the potential and kinetic parts of the symme¬ 
try energy. Assuming that the respective A-meson and 
nucleon-meson coupling constants are the same, the crit¬ 
ical density for the first A - (1232) to appear is found to 
be p c ffff =(2.08 ± 0.02)p o using RMF model parameters 
consistent with current constraints on all seven macro¬ 
scopic parameters characterizing the EoS of ANM at po- 
We also found that the composition and the mass-radius 
relation of neutron stars are significantly affected by the 
formation of A(1232) resonances. In particular, the im¬ 
pacts of the A(1232) formation depend sensitively on the 
values of the g p a and the in-medium Delta mass ttia 
which are also being probed with terrestrial laboratory 
experiments. 

To this end, it is interesting to note that since the 
early work by Kubis and Kutschera [86j], the isovector, 
Lorentz-scalar 5 meson has been incorporated in several 
RMF models in predicting the EOS of dense neutron-rich 
matter and its applications in understanding properties 
of relativistic heavy-ion collisions and neutron stars [87l — 
|95| . Since the <5 contributes negatively to the symmetry 
energy in contrast to the p meson, to fit the saturation 
properties of nuclear matter, reproduce the known sym¬ 
metry energy of about 31 MeV at saturation density and 
avoid the bounding of pure neutron matter, the p-nucleon 
coupling has to be readjusted significantly compared to 
its value in the RMF model without considering the 5 
meson since there is no known mechanism to determine 
the 5-nucleon coupling independently. The net effect of 


including the 5 meson is thus to stiffen the symmetry 
energy at supra-saturation densities where the p contri¬ 
bution dominates. Indeed, applications of these models 
have found that the S may play an appreciable role in 
determining the hadron-quark phase transition [91, 92], 
the crust-core transition and the formation of a pasta 
phase 193, and the mass-radius relation of neutron stars 
[89l l94| as well as the cooling mechanism of protoneu¬ 
tron stars [95|. We have no physical reason to preclude 
possible S meson effects in determining the critical den¬ 
sity for A(1232) formation in neutron stars. However, to 
consider effects of the S meson one has to introduce all 
the necessary but unknown 5-baryon coupling constants. 
Given the fact that our main findings in this work al¬ 
ready depend on the completely unknown A -p coupling, 
the inclusion of the 5 into our current model is thus not 
beneficial for the main purposes of this work. 
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Appendix A: Derivation of equation (12011 

We start from the relation = 2 p n — p p , where 


[—2 j *,2 , — —(3) 

Pn = V (fcp) + ”\iirac + - 5pNPo > 

Pp = \J (^f) I 2 + m dirac + + ffpNPo^, 

PA- =\J (fcp Y + ( to a - g<jA°) 2 + guj AWo - 3g P APo 3) 


(Al) 

(A2) 

(A3) 


are the chemical potential of neutron, proton and A , and m*[ irac = ton — So-Ntr. The condition to determine 
suggests that we can set p&- = 0, kp = 0. Then using the relation above we find that 


m A - gaAcr + goj awo - 3p p APo 3) = \/(fcp) 2 + m^ ac + - g P NPp 3) + Pn ~ Pp, 

where the last two terms on the right hand side can be expressed as [13, [82] 

Pll P P — 4 F/ S y m (/?)5. 

The expression for the symmetry energy in the nonlinear RMF model can be written as [32 ] 

kp 1 


E sym (p) = E^ m (p)+Ef y t(p) = 


9 2 p nP 


6y/k% + m^ac ' 2TO P + Av ^Nff-N^o' 


(A4) 


(A5) 


(A6) 
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/o') 

And the equation of motion for p p -field is given as 

m pPo 3) = SpN (-P<5 - A v PpNPo 3) ^NWo) ■ 
Rearranging this last equation we find 

_(3) = _ g P NpS 

m 2 p + A v g pN g 2 N uj% 


(A7) 


(A8) 


This in turn can be rearranged to give 

-2«?pN^ 3) = 2 jf pNPS 2 -2 ^ 4£^(p)<y = 4(£’ sym (p) - E^ m (p))S. (A9) 

"Ip + Av5pNff,i N Wo 

The wo in the expression of symmetry energy (IA6I) is evaluated at 5 = 0, whereas the one appearing in (IA8I) has the 
value of 5 generally being nonzero, hence we used symbol in the second line. Expanding now [(fcp) 2 + nr^ 2 ac ] 1//2 
as 



*,2 ^ 
dirac 


^ TO dirac 


+ m dirac> 


and using the relation above we find 
/c n, 2 

—^-=$ A + g t tn(1 - x a )a - g^ N (1 - x u )W 0 + p pN (l - x p )pg 3) - 4x p E^ n (p)5 - 4(1 - x p )E sym (p)5 

^ m dirac 

^$A + 5c7n(1 - x a )a - g uN (1 - x u )uj Q + gpvTp^ + ^x p E^(p)5 - 4 E sym (p)6 
+ <?o-n(1 - x a )a - g u N (1 - x u )u 0 + (6x p - 2)Ef°^(p)5 - 4E sym (p)5 
=4>a + 5<tn(1 - Xcr)a - g u n(1 - x^ZJq + (6x p - 6)E sym (p)5 - (6x p - 2)E^ y ^ n (p)5 
=4>a + 5o-n(1 - x a )a - 3 w n(1 - x^ZJq + (6a; p - 6 )E£°£ 1 (p)6 + 4E^(p)5, 


(A10) 


(All) 

(A12) 

(A13) 

(A14) 

(A15) 


where x$ = g^n/g ^a with <fi = cr, lo and p , and $a = m-A — ixin- This last equation is the Eq. EU1) reported in the 
main part of this paper. The derivation of equations (1211) . E^l) and (E3l) follows using similar arguments and steps. 


Appendix B: Derivation of equation (1271) 


Once again we start from the relation pa- = 2p n — p p , and have 

PA~ = 2/ln Pp = Pn T Pn Pp — Pn 4“ 4E/ S y m (p)5 0(5 ) 


(Bl) 


where E sym (p) is the nuclear symmetry energy, <5 = 1 — 2cc p is isospin assymetry, and x p is the proton fraction. The 
chemical potential for neutrons can be expressed as 


_ 9en __ 9£n 2-rf p p d£ N 

dp n dp p 2 85 


where p and <5 are two independent variables. In obtaining the expression above we have used the relation 


(B2) 


8 dp 8 85 8 8 2 r/pj 8 

dpj dpj dp dpj 85 dp p 2 85' 


where J={n, p}, p „ = p p , pp- = p n . When the density is smaller than p c £2, there are only neutrons, protons, electrons 
and muons present in the system. Neglecting the contribution from electrons and muons, £n can be taken as the 
energy density of nuclear matter only. That is, 


(B4) 


e N (p, 5) = [E(p, 5) + m N ]p 
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where ton = 939 MeV is the static mass of the nucleon and E(p,5) is the equation of state of asymmetric nuclear 
matter, which in turn can be written as 


Then, 


E(p, 6) ~ Eo(p) + E sym (p)S 2 + 0(6 4 ). 


(B5) 


—Pn T 4E/ Sym (p)(i 


dE(p,5) , 2t 3 > p d 


—E(p, 6) + m-N + p——~ H-15^" HI ( pE sy m(p )4 2 ) + ^E sym (p)5 


dp 


p 2 36 


— Eo(p ) + E sym (p)S~ + m N + p —^gj ('®sym(p)^ 2 ) + 4-E sym (p)(5 


—Eq (p) + i? sym (/3)5“ + Ton + p 


3 E(p, d) 
dp 


T 4E/ sym (p)$(l £p) 


—Eo(p) + E sym (p)6 2 + ton + p ——■ + E sym (p)(2S + 2 6 2 ) 


=Eo(p) + to n + p 


dp 

3E 0 {p) d , 2 ; 

—~ + T ? {E -" {p)5 > 


+ E sym (p)(26 + 3 6 2 ) 


—Eq(po) + ton + 77 X 2 + 7r~X K 0 + -x 3 + -f-x 2 J( 


3po 


6 ' 


6po 


+ (25 + 3<5 2 ) E sym (po) + ^(26 + 3i5 2 ) x + L(po), 


(B6) 


which is the desired Eq. m■ Note that Kq = Kq(po) and Jq = Jo(po ), and also in the derivation above we have 
made use of the following relation, 


d\ = d_P_Po = J_ 
dp dp 3po 3p 0 ' 


(B7) 
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